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Abstract 

Recently Kaniel & Itin proposed a gravitational model with the 
wave type equation [□ + \(x)]$ a — as vacuum field equation, where 
& a denotes the coframe of spacetime. They found that the viable 
Yilmaz-Rosen metric is an exact solution of the tracefree part of their 
field equation. This model belongs to the teleparallelism class of grav- 
itational gauge theories. Of decisive importance for the evaluation of 
the Kaniel-Itin model is the question whether the variation of the 
coframe commutes with the Hodge star. We find a master formula 
for this commutator and rectify some corresponding mistakes in the 
literature. Then we turn to a detailed discussion of the Kaniel-Itin 
model, file kaniel21.tex, 1998-01-12 
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1 Introduction 

We were very much surprised when we learned during the 8 th Marcel Gross- 
mann Meeting in Jerusalem jjj that Kaniel & Itin || were able to propose a 
gravitational model which looks viable at a first sight even if it had neither 
an Einstein-Hilbert type of Lagrangian nor the Schwarzschild metric as an 
exact solution. Their gravitational potential is represented by a quartet of 
1-forms i? 1 , d 2 1 1} 3 or, for short, by -d a , which constitutes the coframe 
field of spacetime. Their vacuum field equation is simply the wave equation 
with an additional 'massive' contribution depending on some scalar field 
X(x): 

[□ + X(x)] d a = . (1) 
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They show that the Yilmaz-Rosen metric ||, |j solves the tracefree part of 
(|j) exactly. 

Let us be a bit more specific: The Yilmaz-Rosen metric, in isotropic coor- 
dinates, is given by 

g = dt 2 - e 2 ^ (dx 2 + dy 2 + dz 2 ) , (2) 

where r 2 := x 2 + y 2 + z 2 . If we introduce an orthonormal coframe, 

9 = o af3 iT ® d p with o aP = diag(+l, -1,-1,-1) , (3) 

then the following coframe, up to arbitrary local Lorentz transformations, 
represents the Yilmaz-Rosen metric: 

# = e -f- dt, = da; , = dy , and tf £ = dz . (4) 

The tracefree part of (Q) will be determined in Sec. 5.2 and turns out to be 



□ -I( e/3 j D ^) 



tr = o. (5) 



The coframe @j solves the tracefree field equation (||) exactly. We have ver- 
ified this by means of our computer algebra program kaniti . exi displayed 
in the appendix in Sec. [b]. 

Kaniel & Itin tried to derive the field equation (|l|) from a suitable La- 
grangian. For that purpose they had to assume specifically that the varia- 
tion <5t? a of the coframe d a commutes with the Hodge star: *5d a = S*^ 01 . 
However, such a commutativity is only valid for internal Yang-Mills fields. 
It is violated for the coframe and the metric. Therefore the Kaniel-Itin 
model is based on somewhat shaky foundations. 

In the light of the results mentioned so far, the following questions come to 
mind: (i) What is the source on the right hand side of the field equation 
(1)? (ii) Can the Yilmaz-Rosen metric also be adjusted to the trace part of 
(1) and, more generally, to a possible source term on the right hand side of 
(1)? (iii) Is there a consistent variational principle available which would 
allow to derive (|l|), including a source term, from a suitable Lagrangian? 
(iv) What is the (geometrical?) meaning of the constrained variations of 
Kaniel & Itin? 

The purpose of this article is to try to answer these questions. Moreover, 
along our way, we will discuss some unclear points on the commutativity of 
variation and Hodge star which led to some (so far uncorrected) mistakes 
in the literature. — 
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In Sec. g we provide some background material on how to derive wave equa- 
tions of the general type (jl]) from Lagrangians in Maxwell's theory and in 
theories of other internal fields. Here and in the following internal fields 
are those which do not depend on the spacetime geometry (in contrast to 
"& a and ga.fi)- In this way we are able to understand how the Lagrangian of 
Kaniel-Itin comes up in the first place. But since they identify the gravi- 
tational potential with the coframe, we run into trouble from the point of 
view of finding a suitable Lagrangian. 

Any reasonable gauge approach to gravity contains in some way the gauging 
of the translation group. The simplest gauge theories of gravity are telepar- 
allel theories with only the translation group as gauge group. They already 
require the knowledge of how to vary the Hodge dual of forms. In a telepar- 
allel theory, spacetime can be described by an orthonormal coframe i? a as 
the only gravitational field variable, which is interpreted as translational 
gauge potential, see ||. And this gauge potential was used by Kaniel-Itin 
in their model. Accordingly, their Lagrangian is a special teleparallelism 
Lagrangian with the additional postulate of constrained variations. The 
commutativity of 5 and * is, in general, not fulfilled for gauge theories of 
external (or spacetime) groups, i.e., for gravitational gauge theories. In this 
case it is important to know the commutator S* — *8 of the variation 6 and 
the Hodge star *. 

Therefore, in Sec. || we will derive the master formula (^) for 6* — *S. We 
will include general variations of the components g a /3 of the Riemannian 
metric g besides those of a (not necessarily orthonormal) coframe -d a . If 
we insist, in accordance with the Kaniel-Itin postulate, on commutativity 
of 8 and *, then the variations 8g a p of the components of the metric are no 
longer independent and can be expressed in terms of the variation S"d a of 
the coframe, see (|35|). 

In Sec. |] we give a short overview of teleparallelism theories and the rele- 
vant quadratic Lagrangians. We will discuss the viable set of Lagrangians 
and display the results in Table 1. We will show that the KI-Lagrangian, 
for arbitrary variations, is not viable. Some errors in the literature (see 
Schweitzer et al. §, @j) are rectified. 

In Sec. |, we evaluate the model of Kanin & Itin Qj. The field equation 
of the constrained variational principle is the antisymmetric part of a wave 
equation for i9 a , in contrast to the full wave equation as claimed by Kaniel 
and Itin. 

The Yilmaz-Rosen metric, found by Yilmaz Eqs.(18) and (20)] in 1958 as 
a solution in the context of a scalar field theory of gravitation, also turned 
out to be a solution of the bi-metric theory of gravitation of Rosen B; 
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cf. also 



And, in the Kanicl-Itin model, it solves the traccfrcc wave 



equation. In Sec. 5.3 we compare the Yilmaz-Rosen with the Schwarzschild 
metric and give, in Sec. 5.4, a motivation for the emergence of the Yilmaz- 
Rosen metric. Finally, we investigate the implications that would arise if 
the Yilmaz-Rosen metric is considered to be a solution of the field equation 
of Kaniel & Itin including its trace. In Sec. 6 we collect our arguments. 



2 Prolegomena to the Kaniel-Itin model 

2.1 Maxwell's theory and the wave equation 

The kinetic part of the Lagrangian of a Yang-Mills theory is conventionally 
built from the first derivative of the gauge potential A and the corresponding 
Hodge dual. For an internal gauge group, such as for the U(l) or the SU(2), 
the gauge potential A is independent of the metric g or the coframe "& a of 
the underlying spacetime manifold. Then the variation 5 of A commutes 
with the Hodge star operator *. Let us illustrate this for Maxwell's theory, 
i.e., for £/(l)-gauge theory in Minkowski spacetime. 

The Maxwell Lagrangian is given by[| 

^Max = ^ dA A *dA . (6) 

The variation of the 1-form A is independent of the variations Stf" or 5 gap; 
furthermore, it commutes with the exterior derivative, since the variation is 
defined in this way. Therefore, with the coderivative d) := —*d*, we find 

«JL Max = d(6A A *dA) — 5AA *SdA . (7) 

Thus the vacuum field equation reads: 

-*dUA = 0. (8) 

Additionally, we take the Lorentz condition 

SA = (9) 

as a gauge condition. Then, introducing the dAlembertian 

□ := dU + dd* = -*d*d-d*d* , (10) 

1 Wc are using the calculus of exterior differential forms, cf. Our conventions 

are fixed in [Q. 
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the vacuum field equation can be rewritten as 

-D*A = 0. (11) 



One could try to derive ( |11| ) directly by supplementing (|6j) with a suitable 
Lagrangian. The choice 

Ll^:=l^AA^A (12) 

looks suggestive. It leads to the field equation 

-*d^A = 0, (13) 

which should be compared with (||). Consequently the sum of the La- 
grangians (||) and (fl2|), enriched by a matter Lagrangian, 

+ L mat = — (dA A *<L4 + SA A *dt a) + L mat (14) 

would yield directly the wave equation: 

n*A= 6 -^. (is) 

5A v ; 

However, the Lagrangian ( |l2| ) is not gauge-invariant: For the regauging by 
means of the arbitrary function /, 



one findsp 



A — >A + df, (16) 



4a* — 4f« + *dUf A (SA + \d)df} . (17) 



Accordingly, the Lagrangian (14) has to be rejected. We can obtain the 
wave equation ( |i"l| ) only in the special gauge (^|), after the derivation of the 
field equation by means of the variational principle with the Lagrangian (J6j) . 



2 Incidentally, as pointed out by Obukhov, the Lagrangian (jl^) represents an example 
for a Lagrangian of a gauge theory of p-forms, see QQ, Then, instead of @, one has 
A — > A + d* </> as gauge transformation, since 

dtrft = 0/~And the new "Lorentz condition" 

is dA = 0. 
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2.2 A quartet of massive one- form fields 

If we used massive fields, then we would have no difficulties with lack of 
gauge invariance, because the mass term is not gauge invariant anyway. 
Since we want to study gravity a la Kaniel-Itin, we start with a quartet of 
1-form fields k 1 , where I is an internal index with 7 = 0,1, 2, 3. We again 
derive a wave type equation as in the last subsection, but we now add a 
massive term for each of the four fields: 

L k = - (dk 1 A *dk I + Sk 1 A *d)k ! - m {1) k 1 A *k z ) + L mat . (18) 
We vary with respect to k 1 and find as the Euler-Lagrange equation: 

(□ + m (/) )V = ^. (19) 

One could also think of an additional Higgs-type (or 'cosmological') term. 
Then we would have 



L k , = - (dk 1 A *dk' + d^k 1 A *dfik I - m (I) k 1 A *k') 



7f tUKL k 1 A k J A k K A k L + L mat , 



(20) 



and, as field equation, 



(□ + m {1) ) *h' + A eiJKL k J A k K A k L = ^ . (21) 

This is as near as we can approach the field equation (|l|) . Since currents are 
3-forms, we take the Hodge dual of (|l|) and remember * □ = □ *. Further- 
more we put a source term on its right hand side. In gravitational theory, 
this can be only the matter current £„, representing the energy- momentum 
flux of matter. Then the completed Kaniel-Itin field equation reads: 

[□ + A(z)]*tf a = ^=:£ Q . (22) 

Since X(x) is a function, it cannot be identified with some constant mass 
m(7). Also an interpretation of X(x) as a cosmological constant is obviously 
meaningless. Therefore the equations ( pl|) and j22| ) have to be carefully 



distinguished. In future, we will refer to (22) as the (completed) Kaniel-Itin 



field equation. Eq.(E2J) represents the heart of their theory. 
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2.3 Relation to teleparallel theories 

Kaniel-Itin proposed in (|, Eq.(8)] the following Lagrangian for the deriva- 
tion of the vacuum version of ([22]): 

orig 1 

Vki = g A * M <* ~ d)d a A *dH a + A(jc)(i? a A *tf a - 4r/)] . (23) 
We used our notation here. We want to correct this Lagrangian in two 



respects: (i) If we compare (M) with (21), it is clear that we should add the 



first two terms, instead of subtracting them, (ii) The Lagrange multiplier 
term is a hoax since the expression multiplying X(x) vanishes identically: 
We have quite generally d a A np = S^n, with r\p := *i? ( g, and the trace of 
this equation proves our contention. Taking care of both objections, we will 
call 

Vki = \ (d$ a A *dd Q + dH a A *d t $ a ) (24) 

the (corrected) Kaniel-Itin Lagrangian. If we recall that we were able to 
derive (|2l]) from (|o|) only because the k 1 was an internal field, the variation 
of which commutes with the star, then it becomes clear that Kaniel-Itin field 
equation (pl[), for E Q = 0, is not the Euler- Lagrange equation of (fj]|). 

In the Kaniel-Itin model, the Hodge star no longer commutes with the 
variation S'd" since their gravitational potential i9 Q is inseparably connected 
to the spacetime manifold. A postulate of Kaniel-Itin to the opposite, see 
[|[ statement between Eqs.(7) and (8)], is without foundation, at least from 
a geometrical point of view. Consequently, we will give up this postulate. 
The Kaniel-Itin Lagrangian belongs to the so-called teleparallelism models 
of gravity. We will come back to this in Sec. |^. 

The addition of the adjoint piece in ( |2~4| ) docs not break the translational 
mvariance since d*tf a , and thus dH a , can be rewritten in terms of dd a 
(shown here for an orthonormal coframe): 

(TiT =dn a = dd An al3 . (25) 

An analogous procedure is not available in the Maxwellian case for d*A, 
that is, d* A cannot be expressed in terms of dA. Therefore ij^^ is not 
gauge invariant, see (0), and has to be rejected as a decent Lagrangian. 



3 Variation of the Hodge dual of a form 
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3.1 The master formula 

In order to install the variation 6 as a derivation, we demand that it fulfills 
an even Leibniz rule, 

S(u)i A UJ2) — 5lo\ A u!2 + i0\ A 5u>2 , (26) 

where U\ and w 2 are arbitrary exterior differential forms. The Leibniz rule 
is even, because the variation does not change the degree of the form. In 
contrast to this, the interior product v J (here v is a vector) and the exterior 
derivative d decrease or increase, respectively, the degree of the form by one 
and fulfill an odd Leibniz rule. 

Furthermore, we need a relation between the variation and the exterior 
derivative. According to the definition of the variation, they simply com- 
mute: 

[d,S] = 0. (27) 

Let us now turn to the Hodge star operator, see Q| . It maps a p-form 
ip = -if Tp ai ---a p '& ai A- ■ ■Ai3 ai ' into an (n— p)-form ; here n is the dimension 
of the manifold, i.e., in our case n = 4. In terms of components we have 

e ai - ap01 -0 n - v V> 7l .- 7p 0* A • • • A , (28) 

where e is the Levi-Civita symbol. Besides the $-basis |l,^ Ql ,?9 Ql A 
t9 Q2 , . . . A d a2 A ■ ■ • A z9 Q " j, having the Hodge star at our disposal, 



we may define the so-called 77-basis: 

j*l ) *0 otl ,*(i? otl Af»),...,*(r A 'd a2 A • • • A i3 a ' 1 ) j- 



aia2"-a n ^ 

(29) 



Now we can derive the desired expression for 5*0, for an arbitrary p-form 
<f). We can saturate the (n — p)-form *cj) with coframes such as to arrive 
at the n-form 



#fh A • • • A A *4> Q <j> A * (0* A • ■ • A ^) S A rf*"**' . (30) 
We vary (BO). Then the even Leibniz rule (Eq) for the variation leads to 
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S A • • • A A *4> + d 131 A ■ ■ • A A 8*<t> = 

S(f> A rf x "'^ v + cf> A Sr)P 1 '"P p . (31) 

Thus, apparently, we know how to vary the Hodge star, provided we know 
how to vary the ry-basis. The vari ation of the (n — p)-form ?y Q1 "' ap is com- 
puted in the appendix in Sec. A^2. It turns out to be 

(32) 

Incidentally, for the special choice of an orthonormal tetrad, we have Sg a p = 
0, and the last two terms vanish; in particular, we then have SrjPi—P* = 0. 
But we will not introduce this specialization at the present stage. 

If we res olve (31) with respect to 5*<j), then, after some intermediary algebra, 
see Sec. A. 3 of the appendix, we find for arbitrary p-forms </> the master 
formula: 



(6* -*S)$= Sd a A (e a \ *(j>) - * [<5iT A (e Q J 4>) 

1 



Sgap 



t9 (a A(e /3) J» 



■9 



a/3 *, 



(33) 



Again, for orthonormal (co)frames, the terms in the second line vanish since 
5g a p = 0. 



3.2 Constrained variations a la Yang-Mills 
Proposition: The condition 

5*<f> = *S(/) (34) 

for an arbitrary p-form <j> is equivalent to the following relation between the 
variation of the metric and the coframe: 

Sg a p = -2g 7(Q e /3) \Sd' y = -2w( Q/3 ) , with Stf 1 = us 1 V s ■ (35) 

Therefore, for an orthonormal coframe, the allowed variations are of the 
Lorentz type, i.e., U( a 0) = 0- 

To proveQ this equivalence we first assume that (^J) is valid. We apply this 
constrained variation to the volume n-form r\ := *1, defined in (p9|), 

5r) = 5(*1) = *(5l) = , (36) 

3 It was Yuri Obukhov who suggested essential parts of this proof to us. 
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since the constant 1 is not varied. In turn, for the identity 
we find: 

6tf a A np + i? Q A Svjp = . 
The commutation rule ([54|) applied to ftp yields 

§ m = 5(*0 fi ) = *(5tip) . 
Thus, by exterior multiplication with $ Q we arrive at 

0° A Snp — $ a A *{Sdp) = 5dp A *iT = S(gp^) A if . 
On substitution into (|38|), we find 

<W Q A rip + gp-y An a + 8gp 1 An a =0 
or, since np = ep\n, 

$g a p = -2g 7(ct e /3) J(5?9 7 . 
The 1-form 5d a can be expanded with respect to the coframe: 

<5iT = wfti 5 . 
We insert @ into ©. Then we find 

<5ffa/3 = — 2W( a/ g) . 

To investigate the reverse part of the proposition, we apply the gener 
(IH) for the variations of Hodge dual forms and use (|43| ) and (|j) : 

(5* - *S) cj> = uj? a A (e a \ *</>) — * [dp A (e a \ <j>) ] ' 
- 2oj^ d a A (ep\ *<f>) + w 7 7 *(/) 

$a A ( e/3 J*0) — * A (e Q J 0) ] + 5q/3 *<f) 

s^J(^aA*^)-*[^A(e a j0)] 
Thus the proposition is proved. 



/3a 



= 
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4 Teleparallelism theories of gravity 

A Minkowski space is invariant under rigid translations. In order to make a 
manifold locally translation invariant, one can introduce a gauge potential 
by means of a dynamical coframe $ a , see 

spacetime carries a torsion, but no curvature: It is a so-called Weitzenbock 
spacetime. Then, by picking a suitable frame, the connection T a ^ can always 
globally be transformed to zero. 

Therefore, in a teleparallel theory, the coframe d a is the basic gravitational 
field variable. Furthermore, let a metric g be given of Minkowskian signa- 
ture, i. e. of index 3 (number of negative eigenvalues of the metric). For the 
rest of this paper we choose the coframe to be orthonormal, g := o Q( g $ a <S>?? /3 , 
with o a p — diag(+l, —1, —1, —1), and raise and lower the frame indices by 
means of o a p . 

4.1 The Rumpf Lagrangians 

According to Rumpf [^o| , a general quadratic Lagrangian for the coframe d a 
can be expanded in terms of the gauge-invariant translational Lagrangians 
(£ = Planck length, A = cosmological constant = po/2); 



2P 

K=0 



with 



[0] V = itf Q A*z? Q =77, (47a) 

Wy = diT A *d$ a (pure Yang-Mills type) , (47b) 

[2] V = (d$ a A ?? a ^ A * (d&p A ti ) , (47c) 

®V = (d$ a A A * (ffl a A dp) = d& a A -d (i A ( e/3 J*di? Q ) = 2 Wy , 

(47d) 

[4] V = (d<d a Ad 13 ) A* (d0 p A <& a ) . (47e) 
Since <- 3 'V = 2 Wy, we can always put /?3 = 0. 



Variations in tcleparallel theories 



12 



4.2 The irreducible Lagrangians 

Alternatively, the field strength dd a can be decomposed into three pieces 
which transform irreducibly under the Lorentz group: 



(48) 



Here we defined (in parentheses we are mentioning the corresponding names 
of our computer algebra programs) : 

(1) cftT := d& a - {2) dd a - (3) diT (tentor), (49a) 

WcftT := \$ a A (e \dtf p ) (trator) , (49b) 

( 3 )d,r := \e a \ (0? A d&p) (axitor) . (49c) 

o 

In terms of the numbers of components involved, we have the decomposition 
24 = 16©4©4. 

Then we can write 



2P 



1=1 



(50) 



We substitute (|48|) and (|49|) into (|47]) . Then a comparison between ( J46| ) and 
@ yields 

Pi = g ( fl 2 + 2ai) , p 2 = ^ (a 3 - oi ) , p 4 = ~ (ai - a 2 ) , (51) 



Ol = Pi + Pi 



02 = Pi - 2p 4 , 



«3 = Pi + 3p 2 + Pi 



(52) 



and, additionally ciq = po = 2A. These relations were checked by means of 
a computer algebra program, see Sec. O of the appendix. 



4.3 Field equation 

The field equation of a general translation invariant Lagrangian reads Jl5| 

dH a — E a = S a , (53) 

with 

dV _dV_ _ 5L m&t 



(54) 
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In (p4|), the partial derivatives are implicitly defined by means of the varia- 
tion of the Lagrangian: 



dV dV 
5V = 5$ a A — — h S d-d A — — 
dd a ddd a 



(55) 



If we use the abbreviations (p4|), we don't need our master formula for 
the computation of the field equation. Alternatively, one can take the La- 
grangian (^) together with (|47]) and vary the resulting expression by using 
( |33"| ) with 5g a p — 0. This yields the explicit form of the field equation (|53|), 
cf. Kopczyhski [pT|P|: 

-2fT, a = 2pi d*dd a — 2/?2 $a Ad* (dti* 3 A dp) -2 Pi dpAd* (d a A rf^) 
+ pi [e Q J (tW" A *<%) - 2 (e a Jd^) A *c% 
+ p 2 2dd a A * (dd? A dp) 

+ e a \ (diF A d 1 A * (d-d 13 A tip)) - 2 (ejdi^) A dp A * (gW 7 A 7 ) 



P4 



2(% A * (tf Q A d^) 



+ e a J (i? 7 A d??' 3 A * {dd 1 A dp)) - 2 (e a \dd ) A d 1 A * (dtf 7 A fy) 

(56) 

In the first line of this equation we displayed the leading terms containing 
second derivatives of the coframe. In the remaining terms there enter only 
first derivatives. 



4.4 Decomposition of the Kaniel-Itin Lagrangian 

In order to better recognize the structure of the KI-Lagrangian, one can 
decompose it in its irreducible pieces as well as into the Rumpf Lagrangians. 

The first term of the KI-Lagrangian is exactly the Yang-Mills type 
Lagrangian My, cf. with ©. The term d+i? a A *d^d a = -d*d a A *d*$ a is 
left for scrutiny. Using formula ( p5| ) and other rules, see, e.g., |T^| and Sec. 

4 Kopczyriski denoted the Rumpf Lagrangians by K . We have the following transla- 
tion rules: K 1 = My, K 2 = PlV, K 3 = WV. In Q, the second derivatives of the 
coframe are ex actl y the same (for H? = 1) as those in the corresponding three equations 
of Kopczyriski BlL top of p. 503] . 
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A.l for details, we find 



-d*iP A*d*0 Q = -drf A*dn Q = -M A n af) A * (di? 7 A r? a7 ) 

= -iT A ^ A *di?,3 A * (0 a A i? 7 A *di? 7 ) 

= tf Q A ^ A A * [i? a A * (e 7 Jdi? 7 )] 

= ^ A *d^ /3 A ?T A [e Q J (e 7 Jdtf 7 )] 

= *(% A i?' 3 A (e 7 Jd?9 7 ) = * A (e 7 Jd?? 7 )] A dfy 

= -dtf^ A * [e 7 J A di? 7 )] + dtf 7 A *dtf 7 

= -d$ a A A * (iT A dd* 3 ) + dtf a A *dz9 Q = [1 V - [4] V . 

(57) 

Accordingly, the KI-Lagrangian (E4) can be rewritten as 



Vki 



1 



2^ 2 



2 M a A *di? Q - (d0 a A ¥ 3 ) A * (dtf^ A if 
and we can read off the px coefficients as follows: 

pi =1 + 1 = 2, p 2 = 0, p A = -l. 



(58) 



(59) 



By subtracting the adjoint term we would have p\ = p2 = and p4 = 1, 
i.e., we would get the von der Heyde Lagrangian p2[ . 

The coefficients of the decomposition into irreducible pieces, by using (|52j), 
turn out to be 



2 



a-3 



1 



(60) 



(in the von der Heyde case, we have a± = 1, 02 = —2, 03 = 1). Accordingly, 
the KI-Lagrangian can be rewritten in the form 



Via = ~\ 
with the translational "excitation" 

and the coefficients (ISO 



d$ a A H n 



«2 



03 



(61) 
(62) 



The Lagrangian Vki is not locally Lorentz invariant. Rather, a locally 
Lorentz invariant theory results from the following choice of the parameters: 



0,1 = 1 , a 2 = -2 



«3 



(63) 



This represents the teleparallel equivalent of Einstein's general relativity. 
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4.5 Viable Lagrangians 

The form of a general quadratic Lagrangian was displayed in Eqs. ( |50|) and 
( [i"6|) . Various choices of parameters aj or correspond to various telepar- 
allel theories of gravity. We call a specific Lagrangian viable if it leads 
to a theory which fulfills the following conditions: (i) It has the correct 
Newtonian approximation, (ii) It agrees with the first post-Newtonian ap- 
proximation of general relativity, (iii) It has the Schwarzschild metric as 
exact solution in the case of spherical symmetry. 

The question which parameters yield a viable Lagrangian has already been 
discussed in the literature, see for example |2^, [t], [H], [n]. The result is 
that we have viable Lagrangians for a± = 1 , a 2 — — 2 , — arbitrary 
or p4 = 1 , pi = Q , pi — arbitrary. The arbitrary or pi pieces, re- 
spectively, represents the axial square contribution A A * A of the torsion, 
with A := ^ "t}@ A dfip and ^dd a — e a \A. Deviations between viable 
theories due to different axial pieces only show up in fifth order of the post- 
Newtonian approximation |23)| . Therefore, on a phenomenological level, all 
viable teleparallel theories are indistinguishable. In Table 1 we have listed 
some quadratic torsion Lagrangians. 

It follows from the end of the last subsection, see also (|63|), that only the 
teleparallel equivalent of Einstein's general relativity is both viable and 
locally Lorentz invariant. It has not yet been clearly answered if in this 
context local Lorentz invariance is obligatory or merely an aesthetic feature. 





GR|| 1 


vdH H 


viable 


YM 


YMt 


KI§ 


ai 


1 


1 


1 


1 





1 




-2 


-2 


-2 


1 


3 


4 




l 

2 


1 


arb. 


1 





1 


Pi 











1 


1 


2 


P2 


1 

2 





arb. 











Pi 


1 


1 


1 





-1 


-1 



Table 1: This table lists the aj and the pk coefficients for different telepar- 
allel Lagrangians. GR\\ , spelled out in the first column, represents a viable 
gravitational model, the same is true for the von der Heyde case. Obviously, 
the Kaniel-Itin Lagrangian, in the framework of the conventional variational 
procedure, is not viable. We have KL= YM+ YNf and vdH = YM- YNft . 
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4.6 Schweizer-Straumann-Wipf amended 

Schweizer and Straumann [|| and Schweizer, Straumann, and Wipf in- 
vestigated the von der Heyde teleparallelism Lagrangian |2^| , see also ]2jf . 
In particular they showed that, to first post-Newtonian order, the von der 
Heyde theory predicts the same gravitational radiation loss as general rel- 
ativity. However, they assumed in some places the incorrect commutation 
rule 6* — *S. Therefore some equations in these articles must be corrected. 
We stress that the corrections do not influence their overall results, though. 

We present the corrected formulas of the article by Schweizer, Straumann, 
and Wipf in the numbering used there, but in our notation. The cor- 
recting terms are printed in bold. 

The wrong variation first shows up in the explicit expressions of the canon- 
ical energy-momentum tensors derived from the teleparallel version of the 
Hubert-Einstein Lagrangian, e^, and the difference between the von der 
Heyde Lagrangian and the teleparallel version of the Hilbcrt-Einstein La- 
grangian, Ae Q : 



-d [fit 3 A * (<% A <& a )] — dftP A * {d$ a A dp) (2.12), § 

+ \d {Q a A * (d'&P Atfp)} + -d& a A * (dti A dp) 

+ -e a j (A, A d-0 13 ) A * {d-ffi A -dp) 
2 

+-d6n A -dp A e a j* (i? 7 A dti?) 
2 

-ie Q J (d#P A G p ) A * {d&i A i? 7 ) 

-id^ A tip A e a j* (dtf~< A i? 7 ) , 

Ae a = dd a A (dd? A dp) ~^d a Ad* (dd? A dp) (2.13), § 

-ie a J (d&P A tip) A * (dV A i? 7 ) 

~d-d p A i9/3 A e a j* (dtf~< A i? 7 ) . 

Since the additional terms do not influence the antisymmetric part of Ae ai 
we only need to correct the symmetric part: 



(dd a Ad + dd? A d a ) A * {dd~< A 7 ) (3.1), § 
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2 

+~<Mr A A * (dtf 5 A 6? 
2 



Hence we find the corrected field equation as 
e (A-l), 



(dtia + dd A <& a ) A * (diT A # 7 ) (3.3), § 

A di^ Ai? 7 A e a J * (d-d 5 A 



+ idi9^ Ai? 7 A* (dti s A a ) 5f 
2 



In our units, we have for the A-parameter of Schweitzer et al. A = —2p2 = 
| (1 — 03). Like the terms discussed in JtJ, the additional terms are at least 
quadratic in <j) a p (which is the symmetric part of $ Q/ 3 in the expansion 
iT = dx a + $",3 cfe' 9 ), such that the arguments of §, §3] remain unchanged. 

The formulas @, (4.1)] and [Q, (4.6)] have to be corrected similarly as the 
last two equations, replacing d-d 1 by T 7 . Since these explicit expressions are 
not used in the remainder of |?], §4], the conclusions remain valid therein. 

Finally, the quantities A^, B» v , and C» v (formulas @, (5.4b)-(5.4d)]) in 
the expansion of the quadratic and higher-order terms 

eg = ke a s + {A f3a + B f3a + C^} (5.4a), § 

would need corrections. Since these quantities are not used explicitly in Q, 
we did not display the exact expressions here. 



5 Kaniel-Itin examined 
5.1 Lagrangian 

We come back to the Kaniel-Itin Lagrangian (pi|): 

Vki = 7^2 (d$ a A *dti a + d^ a A *d^ a ) . (64) 
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Its Euler-Lagrange equation can be read off from (|5^) by substituting the 
coefficients (|59|): 

-2i 2 S Q = 2 tip A d * A d^) + 4 d*di? a 

+ 2e Q J ((M* 3 A - 4 (eajtft^) A *(% - 2dfy A * (# Q A dtf") 

- e a J [i? 7 Ad/A* (dtf 7 A fy)] + 2 (e a Jdtf") A $ 7 A * (dtf 7 A fy) . 

(65) 

In the first line, we displayed the second derivatives of the gravitational 
potential. As we already saw in Table 1, this field equation is not viable. 

If we use the the constrained variations in (pi]), then we commute 5 and the 
star * and find by simple algebra: 

£ 2 5 (Vja + imat) = 6$ a A (-One +£ 2 Z a )=0. (66) 

Since the variations are constrained, we have to turn to our proposition and 
to use @): 8$ a = up 01 ^, with u [af}) = 0. On substitution in @, the 
field equation turns out to be proportional to the antisymmetric part of the 
wave equation, 

^AQilfl^AEfl, (67) 

rather than to the wave equation itself. Accordingly, the constrained varia- 
tions also lead to a dead end and we have to turn back our attention again 
to the Kl-field equation (p2|). 

5.2 Decomposition of the field equation 

Let us split the full wave equation into its different pieces. For that purpose, 
we most conveniently start from the decomposition of the energy-momentum 
current as a covector-valued 3-form with 16 independent components: 

S Q =X Q +\da A (e 7 j£ 7 ) + ie Q J (<d" A £ 7 ) . (68) 

Here j£ a is its symmetric traceless part, the second term on the right hand 
side its antisymmetric and the last term its trace part, see [ p"2[ eq.(5.1.15)]. 

We come into better known territory if we decompose S Q with respect to 
the 77-basis: 



S Q = T^ a rip 



(69) 
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The T^a's are the components of the energy-momentum tensor. We can 
'saturate' the 3-form T,p by means of the 1-form 



or 



* A T,p) — * (?T A T 1 p 77 7 ) = Tp * (iT A 7? 7 ) 

= T a p*n = T a p**i 

T a p = * (E/3 A i9 Q ) = e a \*Yif3 . 



-T a p (70) 



(71) 

The analog of (|68|) is, of course, the following splitting of the energy- 
momentum tensor: 



T a p — {T( a {3) - - T 7 7 g a p} + T[ af3 ] + - T 7 7 g a p . 



(72) 



Coming back to (22), after some algebra, we find the following decomposi- 
tion: 

Oria - \$ a A (e 7 JDr/ 7 ) - ^e a \ (z9 7 A D Vl ) = I 2 £ a , (73a) 

^ a A tip A (e 7 J D?? 7 ) = l 2 [a A E ffl , (73b) 
i? 7 A Dr? 7 + 4A(x) n = £ 2 tf 7 A E 7 . (73c) 

We can combine the antisymmetric and the symmctric-tracefree part of the 
wave equation in order to get its tracefree part. The fine splitting ( |73| ) 
simplifies to the tracefree and the trace part of the field equation ( |67|) : 



1 



a Va --e a \ (tf 7 AD?7 7 ) =fp a , 
i? 7 A n?7 7 + 4A(ac) = £ 2 i9 7 A E„ 



Sometimes it may be useful to rewrite (74a) by using 

e a \ (i? 7 A Di h ) = -e a \* {* (i? 7 A Di h )} = -* (* (i? 7 A A d a ) 
= -* (tf 7 A Di h ) A *d a = -* (tf 7 A Oft) r) a , 



thereby finding 



or, alternatively: 



□ + 1* A Unp) 



n-\{ep\ne f> ) n a = e 



(74a) 
(74b) 

(75) 

(76) 

(77) 
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5.3 Yilmaz-Rosen and Schwarzschild 
solution compared 

The Yilmaz-Rosen metric and the corresponding orthonormal coframe were 
displayed in (||) to (^). In order to compare the Yilmaz-Rosen metric with 
the Schwarzschild metric, we transform the former one from the isotropic 
coordinates used in (Q) into Schwarzschild coordinates and Taylor expand 
it: 



r 

For the Schwarzschild metric in Schwarzschild coordinates we hnd: 

«« 2m , , 

Soo = 1 - — (exact) , (79a) 

qq „ 2m 4m 2 8m 3 „ /_,4\ . , , 

9ll L = 1 + — + — + — + ° (tit) ■ ( 79b ) 

Their goo components are equal up to second order. The radial components 
gn begin to differ slightly in the second order. Therefore the Yilmaz-Rosen 
solution is consistent with the classical tests of general relativity and can, 
in particular, describe the post-Newtonian perihelion advance correctly, see 
Synge p3, page 296, footnote 1]. It requires further investigations to decide 
whether the two solutions can be observationally distinguished by strong 
gravity effects in close binary pulsar systems. 



5.4 Yilmaz-Rosen metric motivated 

A viable theory of gravitation should be consistent with the local equiva- 
lence principle. Let us consider an electromagnetic wave of frequency to in 
the gravitational field of a point mass to. The frequency shift due to the 
propagation from a point with radial coordinate r to one with r + Ar reads 
(c = G=l): 

Aw . rnAr 

= AU = =- . 80 

uj r z 

According to Mashhoon |25j] , Yilmaz effectively proposed to extend the 
local equivalence principle to a global one. With this idea in mind, we can 
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tentatively integrate (|8I_ 

'' ri du 



UJ 



dU . (81) 



We solve the integrals and find 



= exp[U(n)-U{r )]. (82) 



Using the Newtonian potential explicitly, U = — ^, and taking the limit 
ri — > oo, yields 

Lo(oo) /"m\ , , Ai(oo) / m\ 

— -— = exp — , and thus A , . = exp . (83) 

uj(r) P V r J ' At(r) V r / v ; 

Since At(oo) is not influenced by gravity, one can directly read off 

i9* = dt or g o = . (84) 

Following the pattern of the components of the Schwarzschild metric, we 
now define gu as the inverse of <?oo : 

g = er^dt 2 - e^dr 2 - r 2 (d9 2 + sin 2 6 dip 2 ) . (85) 

However, the Taylor approximation of this metric (which is not the Yilmaz- 
Rosen metric of (||)) reads: 

„ 2m 2m 2 4 m 3 ,„/„„4\ 

5oo = l- — + — -3— + , (86a) 

„ „ 2m 2m 2 4 m 3 ,„/„,4\ , , , 

5n = i + -r + — + 3— + ■ ( 86b ) 

The goo component differs from g^ already to second order. And the devi- 
ation of <7n from gff is doubled in comparison to that of the Yilmaz-Rosen 
gJi \ Therefore, in order to approximate the experimentally well verified 
Schwarzschild metric in an optimal way, we choose the forefactors of ( |85| ) 
as metric components in isotropic coordinates, which eventually leads to 
the Yilmaz-Rosen metric: 

g YR _ dt 2 _ e lf- ^ x 2 + dy 2 + ^ ^ 
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5.5 Yilmaz-Rosen solution and the vacuum field equa- 
tion 

The Yilmaz-Rosen metric (0), keeping in mind (|J) and (Q), fulfills the 
tracefree field equation ( |74a| ) [or, alternatively, ((76|) or ([77])] with vanishing 
source. For a proof compare the corresponding computer algebra program 
in Sec. [§]. As we saw, the Yilmaz-Rosen solution is consistent with the clas- 
sical tests of general relativity. However, its integration constant m cannot 
be directly identified with the source of a spherical body. 



For that reason, we take recourse to the trace equation (74b). In the vacuum 
case we find 

X(x) = ^(^ AD Vl ) . (88) 

The right hand side of this equation can be easily calculated with our com- 
puter algebra program, see also [2j. Therefore the Yilmaz-Rosen metric 
solves the vacuum field equation (22|) provided the 'cosmological' function 
is prescribed as follows: 

^)=-(^) 2 . (89) 

Since such an ad hoc structure looks too implausible to us, we can change 
horses at this moment: We reject the 'cosmological' function of Kaniel-Itin 
and put it to zero. Then we can mimic this function X(x) by means of the 
energy-momentum trace — € 2 T 7 7 /4, cf. ( [74bj ) and (ffl|): 

T\=(^ e -) 2 • (90) 

This energy-momentum trace is plotted in Fig. [j]. Therefore, the Yilmaz- 
Rosen solution fulfills the field equations 

1 /2m 
On a - -e Q J (i? 7 A Di h ) = and * A Dr^) = - ( e~ 



(91) 



Taking an ideal fluid for the description of matter, then, for vanishing pres- 
sure, p = 0, we have T 7 7 = p. Therefore we can understand the above 
computation as a matter distribution (Q) which can be viewed as (proba- 
bly unphysical) star model. The matter of such a star reaches to infinity, 
but it decreases exponentially. We find the maximum of the distribution at 
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r = m /2, compare with Fig. [I]. Most of the star mass is concentrated inside 
the Schwarzschild radius r s — 2m. The volume integral over T 7 7 yields the 
total mass m of the star: 

oo 

J T\ dV = 4tt J T 7 7 r 2 dr 
o 

x\ — — , dr 



where M is the mass in conventional units. 

As a result, we can interpret the constant m of the Yilmaz-Rosen metric 
as the mass of a star, but this mass is distributed in a probably unphysical 
way. 



167rm 2 _!>» 



o 

=^r dx 87rm f „ , 87rm 



£ 2 



Je*d X =^ = M, (92) 



6 Conclusion 

We posed four questions about the Kaniel-Itin model in Sec. 1. We will try 
to answer them in turn: 

(i) We can put the energy-momentum 3-form of matter on the right hand 
side of the Kl-field equation ([!]) , see ( f22|) : 

[n + \tx)] Va = f S ^=:fX a . (93) 

(ii) The Yilmaz-Rosen metric can be accommodated to (^3|) in the follow- 
ing sense: We decompose ( |93"| ) in vacuum into 

-Dn a + ie Q J A Qqp) = , (94a) 

— * (d^ A Oqp) = X(x) . (94b) 

These equations are fulfilled by the Yilmaz-Rosen metric, provided we 
prescribe the 'cosmological' function A(x) in the following way: 

^) = J*(^A^)=-gr?) a . (95) 

If one used a wave equation as field equation, i.e., if in ( p3| ) one put 
A (at) = 0, then one could find the Yilmaz-Rosen solution for the matter 
distribution T 7 7 = (|^ e -T-) 2 . 
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(iii) There doesn't exist a consistent variational principle for arriving at 



(93). Maybe one is able to find one for the tracefree vacuum equation 



(94a) 



(iv) The constraints of Kaniel-Itin on the variations amount to getting 
rid of the independence of the variations of the metric, provided the 
variations of the coframe are prescribed. 

Is the model of Kaniel & Itin viable? Well, it is presently in intensive 
care . . . And it is a beautiful model anyways. 
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A Hodge star and 7/-basis 

A.l Elementary relations for the star etc. 

We now collect some rules for calculations with the Hodge star, where if> 
and </> are forms of the same degree p (see Jl2|): 

**V = (_i)p("-ri+ ind (9) ^ , (96) 

*t/>A0 = *0A^- ( 97 ) 

The index ind(g) of a metric is the number of minus signs if it is in diagonal 
form. Furthermore one has the useful rules 

*(tpA'& a ) . (98a) 

(-l) ind(9) *($ a A*V) , (98b) 

(-lr'^A*^, (98c) 

(_l)&H-l)+ind(ff)^ _ ( 9gd ) 

Sometimes we also need the formula: 



e a \ ip = 
e a \i) = 

*MV>) = 

*(e«JV) = 



i^A(e / ,jV)=pV 



(99) 
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With these rules one can determine the 77-basis, cf. (|2G 



ry Q1 -^ := * (zT 1 A • • ■ A #* p ) = — —r/ ai - a " ap+1 ... an ^ a ^ 1 A • • • A <& a 



(100a) 

1 

(n — p)\ 

— (n — p)! ^l-A.op+i-ain 17 /\ (7 , 

(100b) 

») ai - a " :=*(tf ai A • • • A ) 

= Jldet g^\g a ^ ■ ■ ■ g a ^ e Pl ... 0n = \ e a — . 

(100c) 

Two helpful rules, which connect the different elements of the 77-basis, read 



V ai - ap „ = e„jv ai - ap , (101a) 
v 

§l i A ?/Ql ... Qp = ^(_l)H J («>i^i-« i -iai +I ...« p _ ( 1Qlb ) 



In case of independent variations of the metric components g a p, we need 
the rules 

5 g^ = -g^g^5g y5 , (102) 
8 [det {g^)\ = det (g^) g a ? 6g af3 . (103) 
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A. 2 Variation of the 77-basis 

With these definition and rules, we can compute a general variation of the 
77-basis, involving the fields $ a and g a p: 

frfi-P, = 1 S (if 1 "**/, +1 ,9„^ +1 A • • • A 

(n — £>)! y 



^^^•••%, + ,..^(^A...A^) + ^ 

x <S ( V /| det .9^l5 Ql/31 • ' ' ff QA e ai ... apf3p+1 ... 0n ) A • • • A <>'•• 



(n-p)! 



+ ^^^ 1 ... fe _^^ 4?7 A... ft _ 1 A / 3 i+1 .. A j §gnx 



(100b) 



A7? /3 1 .../3 p ^ +i+ Q ff «A^ 1 .../3 p 



-^(-irv-s'-epj,/ 

1=1 

I „«A /3l.. /3 P _ Ap I / n/t a 



K A 



g Xp e p \ (0* Arf 1 -*')) 5g K x 



(104) 



Hence 



(105) 
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A.3 Deduction of (B3 



We start with (|3l|). We abbreviate ip p := d 131 A • • • A ^ p and get, using the 
variation ( 105 ) of rj^ 1 "'^: 

tp p A 6*4> = 5<f> A Z 1 -^ — dipp A*cf> + (f> A W A rf 1 -^^ 
+ <fi A yd K A ? y ^-^ A - i g^fr-ft.^ s 9tiX 

( S=' © ip p A* (5<P) - Sxp p A *<t> + <t> A A (e^ 1 -^) 

+ (M^A ( 5 a %J^-^) - \^p A ff KA V) &9k\ 



: V P A I * (cty) - - V 3 KA ^aJ - ^ M A (e M J^ p ) A > 
+ (-If ( W A A (e M J V P ) + g X %\ (0 A iT) A *Vp <5 5kA ) 
= V P A (* (Scj>) - 1 V 5 KA <5.9.a) + ^ M A (e M J (0 A > p ) 

- (e M J0) A V P " e M J (V P A » + (-l)^A, A (e M J» ) 
+ ^ p At9 a A*(<M0 k ) <^a 
P ip p A* (8<f>) —ipp A* (60" A (e M J 0) ) + V P A W A (e^J *</>) 
- i> P A - *</> g KX Sg KX + ^ A $ x A g Kp e p \ (» 5g KX . 



(106) 



Since ip p := -d 131 A • • • Atf^" is constructed with p arbitrary $"'s, we conclude 
for an arbitrary p-form </>: 



(6* -*6)cf>= 5$ a A (e Q \ » - * [8$ a A (e a \ 4>) } 



+ Sga/3 



i? (Q A(e /3 )j» - ±g a P *(/> 



(107) 



For the special choice of an orthonormal coframe, we have 5g a p = 0. In 
this case the two last summands vanish: 



{S* -*8)<j>= S$ a A (e Q J» — * Si3 a A (e Q j0) 



(108) 
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B Computer algebra program 

The following Reduce program was written with the help of the Excalc pack- 
age, see [^6|, |27j.[] It verifies (i) the decompositon ( |58|) of the KI-Lagrangian, 
(ii) that the Yilmaz-Rosen metric fulfills the tracefree KI- vacuum field equa- 
tion (74a) and (77), and (iii) the validity of Eq. (89). 

7. file kaniti.exi, 1998-01-11, fwh+fg 7„in "kaniti . exi" ; 



load_package excalc$ 



7. 

7» Basic definitions: 
% 



pform psi=0, r=0, lam=0$ 

f domain psi=psi (x,y,z) , r=r(x,y,z), lam=lam(x,y,z)$ 



psi * d t + a*sin(x)* d x, 

(1/psi) * d x + b*sinh(z)*d y, 

(1/psi) * d y, 

(1/psi) * d z with signature (+1 , -1 , -1 ,-l)$ 

frame e$ 
sgn :=-!$ 



coframe o(t) = 

o(x) = 

o(y) = 

o(z) = 



% 

7» Checking the decompositon (58) of the KI-Lagrangian 
7. 



pform v0rumpf4=4, vlrumpf4=4, v2rumpf4=4, v3rumpf4=4, 
v4rumpf4=4, vki4=4$ 

v0rumpf4 := o(a) " # o(-a)/4$ 

5 This program works properly only with a new patch of Excalc, fixing an earlier bug 
on the hodge dual of scalars. Older versions of Excalc need as additional input the 
following function which should be put in after Excalc has been loaded: 
symbolic procedure dualO u; 
(multpf sq(mkwedge ('wedge . basisf orml ! *) , 
simpexpt list (mk! *sq(absf ! * numr x ./ 
absf!* denr x) (quotient 1 2)))) 
where x = simp!* detm!*; 
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vlrumpf4 := d o(a) ~ # d o(-a)$ 

v2rumpf4 := (d o(-a) ~ o(a)) ~ # (d o(-b) * o(b) )$ 

v3rumpf4 := (d o(a) * o(b)) ~ # (d o(-a) * o(-b))$ 

v4rumpf4 := (d o(-a) ~ o(b)) ~ # (d o(-b) ~ o(a) )$ 

vki4 := rhoO * v0rumpf4 + rhol * vlrumpf4 + rho2 * v2rumpf4 

+ rho3 * v3rumpf4 + rho4 * v4rumpf4; 

rhoO := 0; rhol := 2; rho2 := 0; rho3 := 0; rho4 :=-l; 

diff := vki4 - ( d o(a) " #( d o(-a)) 

+ #(d (# o(a))) ' #(# (d # o(-a))) ); 

diff := diff; 

7. diff has to vanish. Note that our check is not for 
'/, a general coframe. In the following we choose the 
7, Yilmaz-Rosen coframe and set a=0, b=0 : 

a:=0$ 
b:=0$ 

r**2 := (x**2+y**2+z**2) $ 

@(r,x):= x/r; @(r,y):= y/r; @(r,z):= z/r$ 

psi := exp(-m/r)$ 

pform dalembertcof 1 (a)=l , dalemberteta3(a)=3, 
kif eqtrf reea3 (a) =3 , kif eqtrf reeb3 (a) =3$ 

dalembertcof 1(a) := -d(#(d(# o(a) ))) - #(d(#(d o(a) ))); 
dalemberteta3(a) := -d(#(d(#(# o(a))))) - #(d(#(d(#o (a) ) ) ) ) ; 

7. 

7o Checking the vacuum field equation (1) of Kaniel & Itin 
7. 

7. lhs of Eq. (74a) 

kif eqtrf reea3 (a) := dalemberteta3(a) - e(a) _| 

( o(b) " dalemberteta3(-b) )/4; 



7. lhs of Eq.(77) 
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kif eqtrf reeb3 (a) := dalemberteta3(a) - ( e(-b) _| 

(dalembertcof l(b))/4 ) " # o(a); 

% Eq. (88) 

lam := #(o(b) " dalemberteta3(-b) )/4; 
% Eq. (89) 

lam - ( - ((m/r**2)*e**(-m/r))**2) ; 
% Eq. (90) 

energytrace := - 4 * lam / ell**2; 
end$ 
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Figure 1: The prescribed matter distribution £ 2 T 7 7 = (^??e ™) . Most of 
the matter is inside the Schwarzschild radius r s = 2m. 



